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Abstract: We study the stationary problem of a charged 
Dirac particle in (2+l)-dimensions in the presence of a uni- 
form magnetic field B and a singular magnetic tube of flux 
$ = 2-KHi/e. The rotational invariance of this configuration 
implies that the subspaces of definite angular momentum 
1 + 1/2 are invariant under the action of the Hamiltonian 
H . We show that, for K-I>lorK-(<0, the restric- 
tion of H to these subspaces, Hi, is essentially self-adjoint, 
while for < K — I < 1, Hi admits a one-parameter family of 
self-adjoint extensions (SAE). In the later case, the functions 
in the domain of Hi are singular (but square-integrable) at 
the origin, their behavior being dictated by the value of the 
parameter 7 that identifies the SAE. We also determine the 
spectrum of the Hamiltonian as a function of k and 7, as well 
as its closure. 



I. INTRODUCTION 

In Quantum Mechanics, observables are realized in 
terms of self-adjoint operators on a Hilbert space. It is 
for these operators that the spectral theorem holds |L|. 
In particular, the dynamics of a quantum system should 
be given by a unitary group whose generator, the Hamil- 
tonian H (usually a differential operator acting on an 
appropriate space of square integrable functions), must 
be self-adjoint. 

In general, physical considerations lead to a formal ex- 
pression for the Hamiltonian, although they can leave its 
domain of definition not completely specified. Usually, 
one can choose a dense subspace of the Hilbert space on 
which H is well-defined and symmetric, but not neces- 
sarily self-adjoint. 

In these conditions, the question is posed of determin- 
ing if the expression found for H has a unique self-adjoint 
extension in the Hilbert space (i.e., if H is essentially 
self- adjoint), or it admits different self-adjoint extensions 
(SAE) (differing in the physics they describe) and, in this 
case, which one corresponds to the physical system under 
consideration. 

A situation of practical interest in which the Hamil- 
tonian admits nontrivial self-adjoint extensions corre- 
sponds to the movement of charged particles under the 
influence of a Bohm - Aharonov singular magnetic flux 
tube pj , like fermions in the presence of cosmic strings 0] 
or non-relativistic spinless quantum particles interacting 



with a thin solenoid ||4j . In references |3| |6j , this problem 
has been analyzed by means of von Neumann's theory of 
deficiency subspaces JJJ. 

This kind of situations have also been studied as a limit 
of a smeared flux, using a <5-function shell magnetic field 
[ 7,1 or uniform magnetic fields confined to a finite tube 
[ 9 , UJ , and a punctured plane |ll| , |l2|l , which leads to the 
consideration of boundary conditions at a finite radius, 
both spectral and local. It has also been of interest the 
study of charged particle states bounded to flux tubes 

The presence of a <5-like magnetic field has also been 
considered in connection with vacuum polarization ef- 
fects in 18 J, to model the presence of a point-like impu- 
rity in a bidimensional system |l9[ , and more recently to 
describe a non-relativistic electron in the presence of a 
uniform electromagnetic field and a singular vortex, as 
a step toward its application to the quantum Hall ef- 
fect pQ]. This configuration can also be relevant to the 
description of quasiparticles in unconventional supercon- 
ductors |I]@- 

It is the aim of this paper to study the behavior of a 
Dirac electron of mass M and charge e constrained to 
live in a (2+l)-dimensional space, in the presence of a 
constant magnetic field B and a singular magnetic flux 
tube <& = 2irK,/e passing through the origin. In so doing, 
we will use von Neumann's theory of deficiency indices 
to determine the existence of nontrivial self-adjoint ex- 
tensions for the Hamiltonian, a problem that, as far as 
we know, has not yet been solved. 

The rotational symmetry of the problem allows for 
studying the action of the Hamiltonian (a differential op- 
erator H defined on an appropriately restricted set of 
smooth functions) in each invariant subspace character- 
ized by a definite angular momentum I + 1/2. We find 
that the restriction of H to the subspaces with n — I > 1 
or k — I < 0, Hi, is essentially self-adjoint, while for 
< K — I < 1 the operator Hi admits a one-parameter 
family of self-adjoint extensions. In the later case, the 
functions in the extended domain of Hi become singular 
(though square-integrable) at the origin, their behavior 
being dictated by the value of the parameter 7 that iden- 
tifies the SAE. 

Finally, we also determine the spectrum of the Hamil- 
tonian as a function of k and 7. 



II. FORMULATION OF THE PROBLEM 

Let us consider a Dirac particle of mass M and charge 
e in a 2 + 1-dimensional spacetime, in the presence of a 
uniform magnetic field B and a singular magnetic flux 
tube <& = 2irn/e passing through the origin (i.e. the flux 
originated in a magnetic field which is null at each point 
of the plane except at the origin, and whose flux through 
every curve enclosing the origin is finite.) 

The wave function of this particle is a two component 
spinor ip satisfying the Dirac equation (we adopt the fun- 
damental units for which h = 1 = c), 



(ifl- M)i/> = 0, 



(1) 



where the covariant derivative is p = — ie flQ. 

We choose the following expression for the vector po- 
tential leading to the magnetic field under consideration, 



.'Qr k 

A= [— + — 

e er 



(3) 



where 51 = eB/2 has units of squared mass and eg is the 
unit vector orthogonal to the radial direction. 

Accordingly, we get for the Dirac Hamiltonian Hjj = 
v$Tff , where H is the dimensionless differential operator 



m 

-ie ie (-d x - j-de 



H 



-) 

X ' 



\d x -±d e -x 

—m 



-) 

X ' 



(4) 



expressed in polar coordinates (x = vf2r, 6) , with m = 
M/^/Q : the particle mass in units of SI 1 / 2 . 

Since H commutes with the angular momentum oper- 
ator, J = —idg + cr 3 /2, the subspaces spanned by the 
two-component spinors of the form 



1>(x,0) = 



X(x) 



e 

J(l+1)9 



£ L 2 (H 2 ,xdxdO), leZ 



(5) 



are left invariant by the action of H. The restriction of 
H to each subspace characterized by I, Hi, can be cast 
into the form 



1 We choose the following representation of the 7-matrices: 

(2) 



3 1 -22.1 

7 — a , j = — la ,7 = la , 



where the a 1 ', i = 1,2,3, are the Pauli matrices. In a 3- 
dimensional space-time, a non-equivalent representation is ob- 
tained by changing the sign of the matrices, 7^ — > — 7^, but 
this amounts to changing the sign of the parameter M, which 
therefore can be considered to take real values. 



Hi 



(A. 

V dx 



\ dx X 

—m 



(6) 



with a = k — I, when acting on two-component functions 
of the radial coordinate, 



i>(%) 



X(x) 



(7) 



where 4>{x),X{x) € L 2 (R + ,2irxdx). 

In order to ensure that Hi be symmetric and well- 
defined we can restrict its domain to 



V(H l ) = CZ°(R+), 



(8) 



the subspace of functions with compact support away 
from the origin and continuous derivatives of all order, 
which is dense in L 2 (R , 2-wxdx). 

To determine whether Hi so defined is (essentially) 
self-adjoint we must compute its deficiency indices in 
the Hilbert space L 2 (R + , 2irxdx), i.e., the dimensions of 
the characteristic subspaces K± of its adjoint, Hj , corre- 
sponding to eigenvalues ±z, 



n± = dim/C±. 



(9) 



In the following we shall show that Hi admits self- 
adjoint extensions for < K — I < 1, being essentially 
self-adjoint for the other angular momentum subspaces. 



III. SELF-ADJOINT EXTENSIONS 

In order to determine the deficiency indices of the op- 
erator Hi defined in the previous section, we must deter- 
mine the deficiency subspaces K±. 

Let us recall that the domain of Hj , T>(HJ), is the 
' /i(x) 



set of functions f(x) 



e L 2 (R + ,2irxdx), 



for which functions g(x) 
exist, such that 



9 4 X ) ] £L 2 (R+,2iTxdx) 
.92 (x) J 



U,H^) = {g,il>) 



(10) 



for any ip <E T>(H{). The adjoint HJ is defined by g = 

H}f. _ 

Taking into account eq. (M) and the expression for Hi , 
eq. (ph, one can easily see that, away from the origin, the 
first weak derivative of /(x) is locally in L 2 (R + , 27rx dx). 
Therefore, by Sobolev's lemma (see Ref. H), /(x) is ab- 
solutely continuous. This allows for an integration by 
parts in eq. (|l0|), which gives 



9i 

.92 



(s + i+z) 



i(d + l = a_\ 

\ dx x J 



h 
h 

(11) 



In conclusion, HJ acts as a differential operator in 
the same way as Hi in eq. (|6|), but on a larger domain 
2?(-ff/)(D V(H[)), consisting of the subspace of func- 
tions of 1/2 (R) 2ttx dx) which are absolutely continuous 
in R+/{0}. 

In accordance with Appendix |Aj, we must now deter- 
mine the subspaces /C± by looking for linearly indepen- 
dent eigenfunctions of the operator Hj corresponding to 
the eigenvalues ±i, ip/ x y Taking into account eq. (|ll|), it 
is easily seen from 



hHL = ±*v» 



± 

(x) 



(12) 



that the first derivative of ip/^ is absolutely continuous, 



'(*) 



as well as its derivatives of all order. Thus, rjrr,-. G 

C°° P| L2(R + , 2nxdx), and the eigenvalue problem eq. 
(O) reduces to a classical ordinary differential equation. 



Then, eq. (12) leads to the following system of coupled 
differential equations for the components, cf>± and x±, of 
the eigenfunctions ip : 



. dx± 

dx 



1 



X± 



dx 



(» 



(±i-m)^±, 
: (±i + m)x±- 



(13) 
(14) 



Replacing \± from eq. 
other component 



(111) in eq. (p"3|), we get for the 



1 

— < 
x 



& 



2(1 -a) +m 2 + l 



Making the substitution 

</>± = e'^x- a F(x 2 ), 
we obtain Kummer's equation |23J 



= 0. 



(15) 



(16) 



, d 2 F 
d{x 2 ) 2 



(x 2 )+[b-x 2 



dF 
d(x 2 ) 



(x 2 )-aF{x 2 ) =0, (17) 



with a = ?^±1 > and b = 1 - a = 1 - n + I. 

This equation has two linearly independent solutions 
p3| , M{a,b, x 2 ) and U(a,b 1 x 2 ), only the later of which 
leads to <j>± € L2((5, oo),2wxdx), with S > 0. On the 
other hand, the condition 0± e L2((0,<5), 27rcc(ia;) re- 
quires < b < 2 (see Q, pag. 508). 

Moreover, the condition that the second component 
(determined by eq. (|14|)) satisfies \± € ^2(R + , 2irxdx), 
imposes < b < 1. This requires^ that k ^ Z , and 



Notice that if « £ Z, the presence of the singular flux 
through the origin amounts to a shift in the value of the orbital 
angular momentum (as can be seen from eq. (pi)), without any 
further consequence, Hi being essentially self- adjoint. For 
brevity, we will not further consider this case in what follows. 



selects the subspace for which I is the integer part of k, 
re— 1 < I < k, as the only one where nontrivial self-adjoint 
extensions exist. 

Thus, for I ^ [k], Hi is essentially self-adjoint, admit- 
ting a unique self-adjoint extension given by the closure 
of its graph (see Appendix ||) . 

On the other hand, for I = [k] we have found one- 
dimensional subspaces /C±, generated by the solutions of 
eq. (12), ip , given in components by 



*U{a 



l U(a 



1 



4 ' 

X± = 
i 2 + 5 



b = 1 — a; x ), 
— im =F 1 



a;x 



(18) 



(19) 



Therefore, n + = 1 = n._, and H\ K ] admits a one- 
parameter (7) family of (essentially) self-adjoint exten- 
sions W, H? ,, which, as explained in the Appendix 0, 
are in a one-to-one correspondence with the isometries 
IA 1 from /C + onto /C_: 



U^ + 



e^~, 



(20) 



With — 7T < 7 < 7T. 

The functions ipi x \ in the domain of H7-, are of the 
form 



V; = Vo + c(V + + e J7 V"), 



(21) 



where V'o € Cq°(R + ) and c e C, the action of H 7, being 
defined by 

Hfcl> = H [K] i> + ci(tP+~ e*1>-). (22) 

In Appendix |B|, it is shown that the functions in the 
closure of the graph of H[ K ] are continuous and vanishing 
for x — » + . Therefore, the behavior at the origin of the 
functions in the domain of the closure of H7-, , T>{H7 , ) , is 

determined by the behavior of ip = ip + + e lry ip~ , whose 
components satisfy 



(7) = (1 + e* 7 ) 



r(a) 



r(a + 



m 2 + l ' 



0(x a ) 



X 



(7) 



[m(l 



^i\ 



z(l 



0*7 



)] 



r(i 



a 



' m 2 -{-5 > 



r( ! 4 

-1+a 1 /~v^l— a~\ 



(23) 
(24) 



0(x l - a ). 

This allows for the following characterization of the 
boundary conditions the functions ip = I I € 2?(_H7,) 
satisfy: 

0X (7) -X0 (7) ]}=O. (25) 



lim < : 

a:->0+ L 



We will use this con ditio n in the next section to deter- 
mine the spectrum of H? , . 



IV. SPECTRUM OF H 



H 



In this section, making use of the boundary condition 
deduced in eq. (p!5|), we will determine the eigenfunctions 
and eigenvalues of H?-, . So we must solve the eigenvalue 
problem 



ET, 



■y ,/, — 



} ip = \ijj. 



(26) 



Notice that, since H? , is the restriction of W ■> to 

' K- \K.\ 



V(H? K ,) C V(H1 k -,), both operators are realized by the 
same differential operator (given in eq. (|6|), with I re- 
placed by [k]). On the basis of an argument similar to 
the one following eq. ( |12] ) , we conclude that we are look- 
ing for C°° solutions of an ordinary differential equation. 
In terms of the components <j> and Xi we § e t the pair of 
coupled differential equations 



i\ 



a — 1 



+ x X = (A -to) 



i(j)' + i ( — h x) (j> = (A + m)x- 



(27) 
(28) 



Once again, the substitution given in eq. (|lq ) (now 
with < a < 1), leads to Kummer's equation for F(x 2 ), 
eq. (|l7|), with a = — ~ A and b = 1 — a. The requirement 
that <j) an< i X belong to L 2 (R + , 2-Kxdx) selects as the 
unique solution: 



A = e~ x ' 2 X - a U{a = (m 2 - A 2 )/4; b = 1 - a; x 2 ) (29) 

> 
2 — a; x 2 ), 



XA=2(A-m) e -* 2 /V- Q x (30) 



U(a = I + (to 2 -\ 2 )/4;b 
behaving, for x — * + , as 

r(a) 



x 



T(a 



m 2 — A 2 



0(x«), 



_ i(A - m) ^_ 1+a r(l-a 

aA o 



r(i 



m 2 -A 2 



■O(o; 



l-a\ 



(31) 
(32) 



So, the condition expressed in eq. (25) implies 
T(a) 



T{a- 



-x 2 



[to(1 



i(l 



e , 7)] r(i ,i, 



^-»" ^ r(1 m ,-l.J l +'' ' 



r(22 f±5) 

r(a) 



which can also be written as 



T{a 



m 2 + l \ ' 
4 ) 



(33) 



G(A) = (A - to) 

(tan (7/2) - m) 



r(a + m 2 /4-A 2 /4) 

r(l + TO 2 /4-A 2 /4) 

r(a + TO 2 /4+I/4) 



r(m 2 /4 + 5/4) 



Pin)- (34) 



This is a transcendental equation determining the eigen- 
values of H7, . The whole dependence on A is contained 
in G(A), on the l.h.s. This function has simple zeros at 
A = to and A = ±V4 + to 2 + 4n, and simple poles at 
A = ±\/4a + m 2 + An, for n = 0, 1, 2, ... (see Figure 1). 

On the r.h.s. of eq. (|3J), ^(7) is a constant depending 
only on ro, a and the parameter 7 characterizing the self- 
adjoint extension of H\ K y It can take all real values with 
7 ranging from — 7r to 71", being /3(7) > for 70 < 7 < 7r, 
and /3(7) < for — 7r < 7 < 70, where 70 = 2 arctan(TO-). 




FIG. 1. Graphics of G{\) for m = 1 and a = 1/4. The 
horizontal line corresponds to a positive value of /3(7). 



It is evident from Figure 1. that the spectrum of H? , 
does depend on 7. If 70 < 7 < n, the eigenvalues lie 
between a zero of G(A) and the nearest pole on its right: 
For A > to 



y/m 2 +A(N+l) <X N < y/m 2 +A{a + N+l), (35) 
with N = 0, 1, 2, . . . and, for A < to, 



- \/to 2 - 4iV < Ajy < — v/m 2 +4(a- AT-1), (36) 

with TV = -1,-2, -3,... 

For — 7r < 7 < 70 the eigenvalues are bounded on the 
left by a pole and on the right by the nearest zero of 
G(A): For A > to 



\/to 2 + A(a + N - 1) < X N < \/to 2 + AN, (37) 
with A^ = 1,2,3,..., 

- \/to 2 + 4a < A < to, (38) 

and, for A < 0, 

- Vto 2 + A(a - N) < X N < -\/to 2 - AN (39) 

with JV = -1, -2, -3, . . . 

Notice that there is only one level with |A| < 

\Zm 2 + Aa. Moreover, the spectrum of H?, is symmetric 
with respect to the origin only for 7 = 7r and (except for 
the eigenvalue Ao = to) for 7 = 70- 



V. SPECTRUM OF H L FOR L / [k] 



l> [k] 



In this section we complete the description of the 
Hamiltonian spectrum by computing the eigenfunctions 
and eigenvalues of Hi for Z ^ [«] . 

As we saw in Section [II, in the present case Hi is es- 



sentially self-adjoint, admitting a unique self-adjoint ex- 
tension given by the closure of its graph. According to 
Appendix [B], the vectors in T>(H{) are absolutely contin- 
uous functions vanishing at the origin. 

We are looking for solutions of the system given by 
eqs. (p7]-p8j) in this domain. Once again, by an argument 
similar to the one employed in Section III, one can see 
that the eigenvectors belong to C°° f] L 2 (R^ , 27ra dx). 

Following the same steps as in Section [V , one obtains 
the solutions in terms of Kummer's functions. It is con- 
venient to write them in terms of the following pair of 
linearly independent solutions of eq. (p~7|) : 



F x {x 2 ) = M(a;b;x 2 ), 
F 2 (x 2 ) = x 2a M(l + a - 6; 2 - t 7 x 2 ), 



(40) 

(41) 



where a — m _7 A and b = 1 — a, with a = k — I ($■ Z - 
see footnote || -). We will consider the cases I < [k] and 
I > [«] separately. 



l< [k] 

For Z < [k] (a = K — I > 1), only F 2 (x 2 ) leads to 
functions 



t>\ = e 



-* 2 / 2 x«M( — A A " +a-l + a- 1 x 2 ), (42) 



XA 



■I'l 



aM{ 



(m + X) 
m 2 — \ 2 



e - x */2 x -l +a x 



a; 1 + a; x 



m 2 - A 2 + 4a 
4(1 + a) 



2 M(- 



-A 2 



+ a + l;2 + a;x 2 ) 



, (43) 



which are in I/2((0, S > 0), 27ra; dx). Moreover, the condi- 
tion 4>\iX\ € L2HS, 00), 2-rrxdx) requires that M(l + a — 
b; 2 — b; x 2 ) reduces to a polynomial, which occurs only 
when 

Tn — A 
l + a-b = K-lA = -n, (44) 



For I > [k] {a = k — I < 0), only -F^ir 2 ) leads to 
functions 



-* 2 / 2 x- a M( m \ X2 ;l-a;x 2 ), 



X\ 



M{- 



i(m — A) 
2(1 -a) 
m 2 — A 2 



e -x 2 /2 x l- a x 



+ l;2-a;a: 2 ), 



(46) 



(47) 



which are in L2((0, 5 > 0), 2irxdx). Once again, 
the condition <f>\, xx € ^{{5, 00), 2irxdx) requires that 
M(a;b;x 2 ) reduces to a polynomial, which now occurs 
when 



A 2 



(48) 



with n = 0, 1, 2, . . . This time, the eigenvalues are given 
by 



A = ± 2y/m 2 /4: + N, N = 0, 1, 2, 



(49) 



In the present case the eigenfunctions do depend on the 
singular flux, but the eigenvalues are independent of k. 

Finally, notice that in both cases (Z < [k] and I > [k]) 
the eigenfunctions obtained vanish at the origin, thus be- 
longing to the domains T>{H{) of the corresponding op- 
erator. 



APPENDIX A: SELF-ADJOINT EXTENSIONS 
OF UNBOUNDED OPERATORS 

In this Appendix we briefly review the theory of defi- 
ciency indices of von Neumann (for an extended presen- 
tation of the subject, see Ref. fy). We first recall the 
definition of the adjoint of a given linear operator. 

Let A be a linear operator defined on a dense subspace 
T>(A) of a Hilbert space H. The domain of definition 
of the adjoint operator A', T>(A'), is the set of vectors 
■ip G H making the inner product (tp,A(j)) continuous in 
(f> € 1)(A). In virtue of the Riesz - Fischer theorem, for 
any such ijj there exists a unique vector \ G H satisfying 
(ip, A<f>) = (x, 0), V0 S V(A). One defines Aty = \- 

A linear operator A is symmetric if 



with n = 0, 1, 2, . . . So, the eigenvalues are given by 



A = ± 2 v / m 2 /4 + K + N, N = -Z, -Z + 1, -Z + 2, . 



(45) 



Notice that both, the eigenfunctions and eigenvalues de- 
pend on the singular flux k. 



,A<t> 2 ) = (A(t> 1 ,<t> 2 ),\J4> 1 ,(t> 2 e'D(A). 



(Al) 



A linear operator A is self-adjoint if it coincides with its 
adjoint A\ i.e. if V(A^) = V(A) and 



A^(j) = A(j},\f(j)€V(A). 



(A2) 



To establish the conditions a closedf] symmetric oper- 
ator must satisfy to be self-adjoint, a few definitions are 
in order. Let K± = Ker(A' { =p i) be the characteristic 
subspaces of A' corresponding to the ±t eigenvalues re- 
spectively. The deficiency indices of the operator A, n± , 
are defined as the dimensions of the subspaces JC± . 

It is worth recalling that a closed symmetric operator 
is self- adjoint if and only if its deficiency indices are zero 
pf. However, if the deficiency indices are not zero but 
equal the operator admits a family of self-adjoint exten- 
sions whose construction can be carried out by means of 
the following theorem [pj: Let A be a closed symmetric 
operator whose deficiency indices n± are equal; then it 
admits a family of self-adjoint extensions which are in 
a one-to-one correspondence with the unitary maps from 
K+ onto /C_. 

In fact, let IA be such an isometry, then the correspond- 
ing self-adjoint extension Au has domain T>(Au) = {t/j : 
ip = (f> + <f) + +U{<j)+)}, where £ V(A), and <p + e JC+. 
The action of the extension Ajj is given by 



A h 



<t>++U{4> + )) = A{<t>) + i<t> + -iLl{<t> + ). (A3) 



This provides a method for constructing the self- 
adjoint extensions of closed symmetric operators with 
equal deficiency indices by identifying each possible uni- 
tary map from K.+ onto /C_ . 



APPENDIX B: CLOSURE OF H L 

In this appendix we will study the closure Hi of the 
operator in eq. (ra), 



H, 



*(£ + f + ^) 



I A. + i-« _ 

\ dx x 

—m 



(Bl) 



defined on T>(Hi) = Cfi°(R + ), a dense subspace of 
L2(R + , 2irxdx). It will be shown that the functions in 
the domain of definition of Hi are continuous near the 
origin, and vanishing for x — > + . 

In order to obtain T>(H{) we must add to the domain 
of Hi the limit points of the Cauchy sequences in T>(H{) 
whose images by Hi are also Cauchy sequences. 

So, let us consider a Cauchy sequence {VVjJneN with 
ip n € T>(Hi),Wn € N, and such that {J?i^»} n eN is also a 
Cauchy sequence. Therefore, given e > 0, 



\\4>n ~ 4> m \\ 2 < £ 

|| W„ - Mil 2 < e 



(B2) 
(B3) 



for n, m sufficiently large. Making use of eq. (Bl), it is 
easily seen that 



3 Recall that an operator is closed if its graph is a closed 
subset of H x H. Every symmetric operator defined on a 
dense set is closable, i.e., has a closed symmetric extension. 



||^(^-Vm)|| 2 = 

(l^f +P(xM 2 + |xf + <7(,)lx| 2 ) ^irxdx, (B4) 



where we have denoted by (f> and \ respectively the upper 
and lower component of (ip n — ip m ), while the functions 
p(x), q(x) are given by 



p(x) = 


(-+x) +m 2 -2 


(B5) 


q{x) = 


Vl-« V 2 


(B6) 





and are 0(x~ 2 ) for x — > + (since we are taking a ^ Z - 
see footnote || -). It is not hard to see that both p(x) and 
q(x) are positive in the interval [0, 5] for some positive 
S. Only p(x) can change its sign in an interval (x\,X2) 
(depending on a and m), with < 5 < x\ < X2 < 00. 
Notice that the integrand of eq. ( |B4P (obtained through 
an integration by parts) could take negative values only 
in (xi,X2), as a consequence of the term p(x)\(f>(x)\ 2 . 

Moreover, for 5 small enough, we can choose K > 
such that p(x) > K/x 2 . Taking into account eqs. (B2) 
and (B4), for a given e > 0, we can write 



n2 xdx< e, I '—^-dx<e, 
o x 



and 



I \ X '\ 2 xdx<£, f \*Ldx<e, 

Jo JQ x 

if n, m are large enough. Therefore, 

{Vx^' n (x)} and {ip n (x)/Vx} 



(B7) 



(B8) 



(B9) 



are Cauchy sequences in L 2 {0,S) (with respect to the 
usual Lebesgue measure), as well as the sum 

{V£<(a:) + tf»(s)/(2>/£)} - { y^M^)] '}• (BIO) 



Let us call <fr(x) = lim,, 
denote its primitive by 



[y/x ip n (x)}' e L 2 {Q,S), and 



®{y)dy, 



(Bll) 



which is an absolutely continuous function M in (0,(5). 
In particular, \l/(x) is continuous in (0,6). 
On the basis of 

!>/£(*(*) -^«(a0)| = 



My)-(VyMy))'}dy 



< 



My) ~ (VyMv)Y\ 2 dy\ I idy 

in \ Jo 

— > 0, for n — » oo, 



(B12) 



we conclude that the sequence {y/x t/) n (x)} converges uni- 
formly to y/x^(x) in (0, 5), and consequently also in the 
metric of L 2 (0, S), 



lim {\fx^) n {x)} — y/x^(x). 



(Notice that ^(x) is the limit of {ip n (x)} in 
L 2 [(0,5),xdx].) 

In addition, it is straightforward to show that 



tf(ar) 



lim 

n — >oo 



i>n{x) 



(B14) 



Then, we conclude from eqs. (B14) and (Bll) that 



lim {y/xip' n (x)} = y/x^'(x), 



(B15) 



in the metric of £2(0, S). 
Therefore, we can write 



\$>'rxdx < 00, 



1*1 



dx < 



(B16) 



This implies that *'(x) ■ %(x) = 1/2 (*(.t) • *(.x))' € 
Li(0,5). 

On the other hand, the components of ^(x), ^f a (x) 
with a — 1,2, are absolutely contin uous functions in 
(e, <5), for e < 8, by virtue of eq. ( Bll ). In consequence 



J' (KM)' dx = * 2 a (5)-*l(e). (B17) 

In this expression we can take the lim £ ^Q+ , proving that 
the continuous function ^(x) has a well defin ed limit for 
x — > + . Moreover, on account of eq. ( B16 ), this limit 
must be zero. 



As a consequence of the previous results, we conclude 
that the behavior near the origin of the functions in 
V(H?a) is dominated by the functions in K.± (see eqs. 

On the other hand, since the restriction of the Hamil- 
tonian to the subspaces with I ^ [k] is, as already men- 
tioned, essentially self-adjoint, the behavior of the func- 
tions at the origin is dictated by its closure, therefore 
being continuous and satisfying the boundary condition 



lim ip(x) 

a:-t0+ 



0. 



(B18) 
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